Charge current driven by spin dynamics in disordered Rashba spin-orbit system 
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<—j | Pumping of charge current by spin dynamics in the presence of the Rashba spin-orbit interaction 

is theoretically studied. Considering disordered electron, the exchange coupling and spin-orbit 
interactions are treated perturbatively. It is found that dominant current induced by the spin 
00 ' dynamics is interpreted as a consequence of the conversion from spin current via the inverse spin 

, Hall effect. We also found that the current has an additional component from a fictitious conservative 

field. Results are applied to the case of moving domain wall. 

Recent spintronics studies aim at manipulation both of charge and spin degrees of freedom Central roles 

■ are played by the spin-orbit interaction and the exchange interaction between the conduction electrons and local 
r/3 \ spins [2, 0, 0, 0, 0] . It has been shown that the exchange coupling is useful for electrical control of magnetization 

dynamics via spin transfer torque 0, 0]- It can a ls° be used to pump spin current from the precession of the 
magnetization 0, 0, 0] . The spin-orbit interaction has been recently found to induce a magnetism by the application 
of electric voltage (the spin Hall effect) [l(J El, E2] • 

By combining the exchange and the spin-orbit interactions, various phenomena are expected, and the subject 
discussed in this paper is one of them; pumping of charge current by dynamical magnetization. The idea is to convert 
the pumped spin current into a charge current by using of the spin-orbit interaction as proposed by Saitoh et al. [3]. 
c"j ' This current due to the inverse spin Hall effect was indeed observed , [HI, EH in metallic systems where the spin-orbit 
O ! interaction is induced by Pt atom. 

Theoretically, generation of the electric field or voltage due to the dynamical spin structure was discussed by 
Stern [l6[. The field is not a real electric field, but an effective one due to a spin Berry phase acting only on charge 
C*") ■ degrees of freedom with spin (like the electron) . The mechanism is similar to the Faraday's law, but a magnetic flux 
is replaced by a fictitious field from the spin Berry phase. The induced field is described by V x E e g = —b, where b is 
field of spin Berry phase, and current is divcrgcncelcss; V • j = 0. The theory was recently applied to a domain wall 
\ by Barnes and Maekawa [l7j], and the effect of the spin relaxation was studied by Duine [lij], where the relaxation 
C^) • was introduced by a phenomenological term (/3-term). All these studies have been done in the adiabatic limit, where 
\ the exchange coupling between the local spin and the conduction electrons is strong. 

■ Charge current as a result of inverse spin Hall effect was theoretically studied by Zhang and Niu [l9| and Hankiewicz 
f*»» \ et al. [20( (the effect was called reciprocal spin Hall effect). The calculation was done as a response to applied spin- 

■ dependent chemical potential, which would not be easy to control experimentally. In contrast, our study tries to 
derive direct relation between physically accessible quantities, current and magnetization (local spin). 

Another type of a voltage generated at an interface of a ferromagnet-nonmagnct contact was predicted by Wang et 
al. 21 1 in order to explain the experimental observation [22| . They have pointed out that a net charging is due to a 

■ back flow of spin current at the interface, and the spin accumulation at the interface is essential. 
In this letter, we theoretically predict another mechanism of a spin-induced charge battery realized in disordered 

conductors. We consider the perturbative regime of the exchange coupling, that is in the opposite limit of adiabatic 
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cases [lg, [17|, [18| . Proposed mechanism does not rely on the interface, and the pumped current arises simply when 
the exchange interaction and the spin-orbit interaction exist simultaneously. 

We consider the two-dimensional electrongas system (2DEGs) with the Rashba spin-orbit interaction that originates 
from the lack of the inversion symmetry [H, 0, [23|] • The conduction electrons couples to the local magnetic moment via 
the exchange coupling. Such a system can be achieved by the 2DEGs attached to the ferromagnetic contact [24, 25|, 2(|, 



or the 2DEGs in magnetic semiconductors (e.g., CdTe/CdMnTe) [27|,[28| 



Let us consider the current representation by the simple argument. The current j^, proportional to the average 
tr(fc^) of the electron wave vector with respect to electron states (tr is trace over spin indices), vanishes if the 
system is spatially symmetric. The Rashba interaction, proportional to (k x <x) 2 , breaks the spatial symmetry, but 



FIG. 1: Diagrammatic representation of currents at the first order in J ox and the lowest order in £7, which turns out to vanish. 
Dotted lines and wavy lines denote local spins S and the Rashba interaction, respectively, and thick line represents the diffusion 
ladder, D q . 



charge current docs not arise since tr(fc /J (A; x <x) 2 ) = (only spin current arises Q). Charge current appears when 
we introduce the exchange coupling, proportional to S ■ it where S is a local spin. Wc would have the current 
j M oc tr(fc p (fe x &) Z S • & ) oc t^vzS 1 ' to the first order exchange interaction, and oc (S x S') v at the second oder 
interaction with different spins, S and S' . 

In order to obtain the precise expression of the current, we perform analytical calculations by using a diagrammatic 
technique. Although the argument in the previous paragraph gives the qualitative idea of pumping charge current, 
it turns out that the linear term in S vanishes identically in the Rashba case, and the second order term needs a 
dynamical part as S x S. We will demonstrate that the pumped current has two components. One describes the 



inverse spin Hall effect [13l ] , and the other is a conservative current which written as a divergence of a scalar potential. 
We will show that the current due to the inverse spin Hall effect is dominant in various cases, and the conservative 
current is relatively small. However, the conservative one is also important from the view of the fundamental physics, 
because it provides the fictitious scalar potential which acts only on the particle having both charge and spin. 

We consider a disordered 2DEGs with the Rashba spin-orbit interaction. The 2DEGs also interact with the local 
spin, S x (t), via the exchange coupling. The local spin is treated as classical, and slowly varying in space and time. 
The system is represented by a Hamiltonian H = Hq + H cx (t) + H so + Hi mp , where Hq = J^ fc Ek,c^ k Ck describes free 
electrons with Ek = h 2 k 2 /2m (m being the effective mass). The exchange and the Rashba spin-orbit interactions are 
given by 

H m (t) = -J cx ]T S q (n)e int cl +q &c k , 

h,q,Q 

H so = -a^e A1I/z fc A1 (4ir iy c fc ), (1) 

k 

where J ox is the exchange coupling constant, S q (Cl) denotes the Fourier transform of the local spin structure 
and a is the strength of the Rashba spin-orbit interaction. Spin-independent disorder is represented by H lmv = 

Sr=i Sfe k' Y e ^ k ~ k ^ *" ic L Cfc > wn i cn gives rise an elastic electron lifetime r = (2-Kpn\v? /V) -1 , where p is the density 
of states, Hi is the number of the impurities, u is the strength of the impurity scattering and V is the volume of the 
system. 



The charge current density of this system is given by 



i.OM) = ^e^-tr 

k,k' 



(fc + fe')„ 

2m -<*W^)G£ fe ,(M) 



(2) 



G^ k ,(t,t') is a lesser Green function which is a 2 x 2 matrix in spin space with components G^ a k , a ,(t,t') = 

i( c Lo-'(^') c fcer W) l* 7 ' G ' = = ' = )' where (■ ■ ■) is the expectation value estimated by the total Hamiltonian H. 

We calculate the current by trea ting both Rashba (to the first order) and exchange interactions perturbatively, which 
is valid if J cx t 1 and afcpr <C 1 [29(. Successive impurity scatterings are denoted by ladder approximation, resulting 
in a diffusion propagator at small momentum transfer (q), D q = (Dq 2 ^^ 1 , where D = k^T /2m 2 is a diffusion constant. 
Dominant contributions are from diagrams that include a maximal number of diffusion propagators. Contributions 
from the first order in J ox are shown in Fig. [TJ that turn out to vanish identically. The leading contribution coming 
from the second order in J ex are shown in Fig. [2l 



FIG. 2: Leading contribution from the second order in J cx . Contributions from (a) and (b) are denoted by and j^ SH in 
Eq. (|4]l. respectively. Contributions from other diagrams cancel out or are smaller by 0(1/eft). 



After straightforward calculations, the current in the slowly varying limit (Sir <C 1) is obtained as 



4eaJ cx 



-e^(^)£e- iQ -*D g (S Q - g xS 9 ), 



(^f)D q D Q A^ q C qQ + D q E$ Q B" q DqAqFqq 



(3) 



where A£ = E fe *^_Q ^ = £fc M.9fc) 2 fffc +<r ^ qQ = Re £ fc ^ Q = 

ilm£ fc (fc + §)^5fc.9l +q .9fc +Q and F v qQ = RcJ2 k {k + q) v g r k {g% +q ) 2 gl +Q (g T k = (g%)* = (e F - £fc + (e F being 

the Fermi energy)). The second and third terms in Eq. (J3j) (proportional to E q Q and F q Q, respectively) are the second 
leading term with a long range limit (q, Q — > 0). It is, however, physically the most essential term as we see below. 
For a spatially smooth structure of spins, i.e., q,Q <C I (1 is the electron mean free path), we can approximate 



jet e* 1 , 

obtained as |3fj| j^{x, t) = j]f K (x, t) + jt(x, t), where 



2nipTmDQ fl , B v q ~ 4:irpT 3 £ F q v , C qQ ~ -/o/2ef, x^ qQ 



-■K'xp^q^ and F£ q 



npr Q v . Then, the current is 
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Dx — xi ( ^ a; ^ ^a: 1 ) r; : 
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D, 



dD\ 



(4) 



V 



(2) 



E g e lqx (D q ) 2 and a is the lattice constant. (Note that singular behavior at 



Here, £>, 

q —> is cut off at g ~ L , where L is a system size.) 

The first term in Eq. (j4j, j ISH , describes a current whose direction is correlated with the magnetization direction, 
perpendicular to S x S, and represents inverse spin Hall effect fl3j | . In order to make clear the physical meaning of 
this current, we compare with the spin currents pumped by magnetization. In the absence of spin-orbit interaction, 
we obtain the pumped spin current at the lowest order in J ex as 



Js M (a:,f) 



Je^FT 2 d f d 2 Xi 



2ir dx. 



o F)x — X\ 



2Jc: 



d 2 x 2 



2 Dx 1 -x 2 



(S Xl x Sx 2 ) 



(5) 



The polarization of the spin current is in both directions, oc (S) and (S x S), where (• ■ •) denotes average over diffusive 
electron motion. This spin current is a gradient of a certain spin potential, j S[i = — V^0 S . The result of Eq. ([5]) is 
consistent with the observation by Tserkovnyak et al. Q, where the spin current appears at the interface between 
ferromagnet and normal metals associated with the phenomenological parameter of spin-mixing conductance. (We note 



that there is also an equilibrium component of spin current [31|. This component, j s < ^ q \x) 



free from diffusion poles. Hence, it is local and therefore small compared with dynamical contributions.) The meaning 
of pumped spin current is understood by taking a divergence: 



V-j B (x,t) = - 



2ir 



-2J n 



cPxi 



Dx — xi ( & x £>x i , 



(6) 
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FIG. 3: (Color online) Two typical geometries with ferromagnets attached to two-dimensional electron system. Left: Datta- 
Das geometry [24|j and Right: perpendicular geometry like in Ref. [l3||. In both cases, inverse spin Hall current is given by 
jj? H ace^S x S)„ (n = x,y). 



Comparing the second term (V ■ j^") to j ISH , we see that j^ SH = 7ish£^i/z(V • jj" 2 ^), where 7ish = —Sear. This 
expression is the Rashba- version of inverse spin Hall effect, j oc j s x <r, proposed in Ref. (Note that the spin 

current considered in Ref. (l3j is the one flowing through the interface that enables the spin current to enter without 
divergence.) Eq. ((6]) represents a conservation law of spin, and correctly describes a fact that the Gilbert- type damping 
((S x S)) results in a flow of spin current or S. 

In contrast, the second term in Eq. ([4J, j^, is a gradient of a scalar quantity. It can be interpreted as a current 
arising from a potential or a conserved force. The fictitious electric field, defined by Ei n d = j^/ao, where ao = e 2 nT/m 
is Boltzmann conductivity, is written as .Eind = — V<^ind- The scalar potential is obtained as 

2aJ ex V [tfx! [d 2 x 2 dD%) ■ 



iree-p J a J a ox\ 

(Note that these scalar potential and field are fictitious ones, acting only on charge having spin degrees of freedom.) 
The current is in the direction where magnetization changes. It contributes to the perpendicular current in the 



Datta-Das spin transistor geometry 24|. However, it is not in-plane current in the layer geometry [13| as shown in 
Fig. [31 It does not contribute to the case of moving domain walls as we will show below. 

Let us apply our results to a case of a moving domain wall (as would be realized by using magnetic semiconduc- 
tors [U Hal)- We define polar angles (9,4>) with respect to the easy and hard axis of spin as cos 9 = S easy /S,ta,n(f) = 
Shard/ S m id, where S ea sy, Shard and S'mid denote spin components in easy-, hard- and medium- anisotropy directions. 
Rigid and one-dimensional (in the x-direction) domain wall solution is represented by two dynamical variables, X(t) 
and 4>(t) [13], as cos9(x,t) = tanh x and sin(?(x,£) = [cosh — y^ ] -1 , where A is wall thickness. By assuming 
a dirty case A 3> £ (I being the mean free path), and noting that S x x S x ' vanishes if \x — x'\ A, we can approx- 
imate S x x S x < by local value as S x x S x ~ sin9(x, t)(<pe s — -f-e^). Here, e e = (cos #cos <j>, cos 9 sin^, — sin#) and 
= (— sin </>, cos </>, 0), indicating that the damping (S x S) on the translational motion (X) is in the (^-direction, 
while it is within the wall plane when <p varies. The pumped inverse spin Hall current, j^ SH ~ — Jo e /jjjz(£ x S) v , where 
jo = ^emaJ cx 2 T 2 Dq (Dq = \ dxD x ) depends much on the wall geometry. We consider three types of the domain 
wall as shown in Fig.HJ a Neel wall (N) and two Bloch walls ((B-i) and (B-ii)). By estimating (S x S) inside the wall 
(by using (cos 9) = (tanh §) = etc.), we obtain 

. ISH (N) . X .ish( b -0 - X . 

3 X jo-^sm0, j£ a -jo— cos0, (8) 



which are driven by translational motion, and 



which is driven by tilt of the wall. In the case of the Neel and the out-of-plane Bloch wall, the current is a constant 
(if wall velocity is constant) at small speed and shows an oscillation (sin</> or cos<^>) when the domain wall is above 
Walker's break down. In contrast, no current is induced for the in-plane Bloch wall at small velocity and finite but 
steady current arises above breakdown (as long as <p is more or less constant). The gradient part of the current, , 
on the other hand, is averaged out for any type of the domain wall. 

Let us briefly see the magnitude of current (Eqs. ©(O). We use £f = lOmcV, J cx = lOmcV, a = lOnm and 
a = 0.3 x 10~ n eVm Q. Using Dq — (L 3 /a 2 £ 2 X) ln(L/A) (our diffusive result depends much on sample size L), the 
current is estimated as j — 4 x 10 _14 [C/m] x In j x fi[Hz]. If we choose fl = 100MHz, L = 1/xm and A — 10a, we 
obtain j — 10[A/m], i.e., current is I = jL — 10/iA, which would be detectable experimentally. 
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Neel (easy axis // x) Bloch-i (easy axis Hi) Bloch-ii (easy axis // y) 
FIG. 4: (Color online) Three domain walls in the xj/-plane, Neel, Bloch-i and Bloch-ii. 



Wc have considered a case of an uniform Rashba interaction. If we allow a to be position dependent. a(x), we have 
other contributions. For instance, charge current linear in S arises proportional to ((Va)S>). Thus, various currents 
are expected in the case of finite-size Rashba system (finite size is always the case in experiments). 

In conclusion, wc have theoretically shown that charge current is pumped by magnetization dynamics in the presence 
of the Rashba spin-orbit interaction. The dominant part was found to be due to the inverse spin Hall effect, i.e., 
conversion of spin current into charge current by spin-orbit interaction. In addition to the inverse spin Hall current, 
we found a conservative current flowing basically along the gradient of the magnetization damping. This current 
is rotation free, and should be distinguished from the inverse spin Hall current and from the divergenceless current 
predicted by Stern and others 0, 17 , [l8| . It would be extremely interesting if one could experimentally determine 



the type of the pumped current. 

The authors are grateful to E. Saitoh, B. Kramer, R. Raimondi, M. Yamamoto, S. Kettemann, J. Shibata, H. 
Kohno, S. Murakami, and T. Ohtsuki for valuable discussions. This work has been supported by the Deutsche 
Forschungsgcmcinschaft via SFBs 508 and 668 of the Universitat Hamburg. 

Noted added in proof. — After finishing the manuscript, we found optically induced inverse spin Hall effect was 
observed in GaAs 13311. 
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